Metal-insulator transition in quarter-filled Hubbard model on triangular lattice and 

its implication for the physics of NaQ^Co02 
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The metal-insulator transition of the quarter-filled Hubbard model on triangular lattice is studied 
at the mean field level. We find a quasi-one dimensional metallic state with a coUinear magnetic 
order competes closely with an insulating state with a non-coplanar magnetic order for both signs 
of the hopping integral t. In the strong correlation regime(C//|t| ^ 1), it is found that the metal- 
insulator transition of the system occurs in a two-step manner. The quasi-one dimensional metallic 
state with coUinear magnetic order is found to be stable in an intermediate temperature region 
between the paramagnetic metallic phase and the non-coplanar insulating phase. Possible relevance 
of these results to the physics of metal-insulator transition in Nao.5Co02 is discussed. 

PACS numbers: 



NaxCo02 is a remarkable transition metal oxide sys- 
tem in which strong electronic correlation and geomet- 
ric frustration interplay with each other. Unconven- 
tional superconductivity, novel magnetic and charge or- 
dered states and state with peculiar Curie- Weiss metal 
behavior are observed in this system at different carrier 
concentrations [ll-Q. To date, most of these observations 
remains not well understood. 

The quarter-filled system, namely NaQ,^Co02 is es- 
pecially interesting!^, Q. It divides the phase dia- 
gram of the NaxCo02 system into two parts of quali- 
tatively different nature. For x < 0.5, the system ex- 
hibits conventional metallic behavior in the normal state 
and with water intercalation a superconducting state 
with unknown paring symmetry in a small doping region 
around x = 0.3[13. While fore x > 0.5, the system ex- 
hibits Curie- Weiss metal behavior at high temperature 
and a weak static magnetic order at low temperature 
with ferromagnetic intra-layer spin correlation 0, d, 
Right at a; = 0.5, the observations are even more strik- 
ing. The system undergoes two phase transitions at 88K 
and 53K[3|, Q. At the 88K transition, the system de- 
velops a magnetic order with a ordering wave vector 
q= i(bi +b2)(9|. However, the longitudinal resistivity 
of the system seems to be totally ignorant of the transi- 
tion. While at the 53K transition, the system enters the 
true insulating state as signaled by the rapid increase of 
the longitudinal resistivity below this temperature. How- 
ever, both the peak position and the intensity of the mag- 
netic Bragg peak developed at 88K seems to be ignorant 
of the transition at 53K. The observation of two transi- 
tions and the apparent mutual independence of the mag- 
netic order and the charge transport have arose many 
interest in the field. However, this is still not the whole 
story. In fact, it is found that Hall resistivity begins to 
drop abruptly around the 88K transition which indicates 
that it is not a transition in the spin sector onlyQ. 

No coherent picture exists yet on the low tempera- 
ture phase transitions of the NaQ c,Co02 system. It is 



generally anticipated that the ordering of the Na ions 
outside the C0O2 plane may plays an important role in 
the low temperature physics of this system 0, @, [ll], 13 1. 
Charge ordering(induced by the Sodium order) in the 
C0O2 plane is invoked in many works to explain the 
multiple phase transitions. In this paper, we suggest a 
possible picture for the multiple phase transitions of the 
Naa,5Co02 system as an intrinsic property of a quarter- 
filled strongly correlated system on the triangular lattice. 
We backup our reasoning with a mean field study of the 
metal-insulator transition of the quarter-filled Hubbard 
model on triangular lattice. We find this seemingly over- 
simplified model does exhibit two successive transitions 
at quarter filling in its way toward the low temperature 
insulating state. In our picture, the magnetic transition 
at 88K transform the system into an quasi-one dimen- 
sional metallic state through a spin ordering induced di- 
mensional reduction process. This quasi-one dimensional 
metal can then transform into the true insulating state at 
53K in various possible ways. We thus predict that the 
intermediate phase between 53K and 88K to be a state 
with one dimensional transport property. 



FIG. 1: (a)The magnetic structure of the proposed state for 
the intermediate temperature region between 53K and 88K. 
Note the state has no charge order and all lattice sites are still 
uniformly quarter-filled. (b)The insulating state with a four 
sublattice and non-coplanar magnetic structure predicted by 
the mean field theory of the quarter-filled Hubbard model 
on triangular lattice. The magnetic moments on the four 
sublattices point to the four corners of a perfect tetrahedron 
as shown in the inset. 

The state we propose for the intermediate temperature 
region between 53K and 88K is shown schematically in 
Figure la. In this state, the system self-organizes into a 
series of chains with ferromagnetic order on them. Neigh- 
boring chains are antiferromagnetically correlated with 
each other. Unlike other proposals, there is no charge or- 
dering in this state and all lattice sites are still quarter- 
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filled uniformly and thus the ferromagnetic chains are 
metallic. It is interesting to note that this metallicity is of 
one dimensional nature as any inter-chain electron hop- 
ping must cross the spin gap. Thus the spin ordering in 
such a particular pattern renders the system a one dimen- 
sional metal and so we dub it as a spin ordering induced 
dimensional reduction process. Similar dimensional re- 
duction process also occurs in other condensed matter 
system with multiple low energy degree of freedoms. The 
most well known example is the orbital ordering induced 
dimensional reduction in ruthenium oxides [l5|. 

The magnetic structure of the state shown in Fig. la is 
given by < Si >= mexp(iQ3-Ri) with Q3 = i(bi+b2). 
The states with Qi = \hi and Q2 = |b2 are degener- 
ate with it by lattice symmetry. All the three states 
are consistent with the neutron scattering measurement 
and the real system will contain domains of all three 
kinds. Since the ferromagnetic chains are still metallic, 
the longitudinal resistivity of the system is not expected 
to change dramatically across the 88K transition. How- 
ever, the transverse transport is totally different. As the 
system becomess quasi-one dimensional in charge trans- 
port state, we expect the Hall resistivity to drop abruptly 
at the 88K transition, in accordance with the experimen- 
tal observation. Thus the state shown in Fig. la has the 
correct characteristics to account for the phenomenology 
of the intermediate phase between 53K and 88K. 

With further lowering of temperature, the quasi-one 
dimensional intermediate state will eventually become a 
true insulator. It is important to note that the inter- 
mediate state has an enlarged unit cell with two lattice 
sites. Thus the system is half-filled below the 88K tran- 
sition. To arrive at the true insulating state, the sys- 
tem should lower its symmetry further to have a unit 
cell with at least four sites. Through mean field cal- 
culation on the Hubbard model on triangular lattice, 
we find a state with a four sublattice structure and 
non-coplanar spin configuration (see Fig. lb) that meet 
this requirement. This state has no charge order either 
and the magnetic structure can be written as < S^ >= 
nil exp(iQi • R^) + xa.2 exp(iQ2 • R^) -I- ms exp(iQ3 • R,). 
Here the three vectors mi , rh.2 , ma are of same length 
and are mutually orthogonal. The magnetic moments in 
the four sublattices points to the four corners of a per- 
fect tetrahedron. This state has the same set of magnetic 
Bragg peaks as the collinear state shown in Fig. la if the 
three kinds of domains of the collinear state appear with 
the same probability. The mean field theory presented 
in this paper also indicates that the transition between 
collinear state and the four sublattice insulating state 
leaves the intensity of the magnetic Bragg peak almost 
intact. This is in accordance with the observation at the 
53K transition. 

It should be noted that both kinds of order shown in 
Fig.l have already been predicted in a mean field study 



of the quarter-filled Hubbard model on triangular lattice 
with a negative hopping integral [l3|. In that work, both 
kinds of order develop as a result of the nesting property 
of the Fermi surface and exist for any non-zero interac- 
tion. However, the Fermi surface of the NaQ,^Co02 is al- 
most featureless. As will be shown below, the Nao,^Co02 
system should be described as a quarter-filled system 
with a positive hopping integral. Since the triangular 
lattice has no particle-hole symmetry, any magnetic or- 
der can develop only for interaction greater than certain 
critical value. Indeed, we find the magnetic structures 
shown in Fig. 1 exist only in the strong correlation regime 
for Naa,^Co02- 

We now present the results of the mean field study of 
the metal-insulator transition in the quarter-filled Hub- 
bard model on triangular lattice. The Hamiltonian of the 
model reads 

H = -t ^ {cl„CJ,^ + h.c.) + uY,n^■:\n^A■ (1) 

<i,j>,<J i 

It is generally believed that the NaxCo02 system re- 
sides in the strongly correlated regime with ^- » l[l3]. 
Since the triangular lattice has no particle-hole symme- 
try, the sign of the hopping integral t is important. For 
NaxCo02 system, it is negative as inferred from ARPES 
measurement fl3| . In Nao,5Co02, each Co site hosts 1.5 
electrons on average is thus a |-filled system with a neg- 
ative hopping integral. After a particle-hole transforma- 
tion, it is equivalent to a quarter-filled system with a 
positive hopping integral. For purpose of comparison, 
we also consider the case of quarter-filled system with a 
negative hopping integral which is studied already in Ref 
12. 

From the weak coupling point of view, a quarter-filled 
system on triangular lattice with a positive hopping in- 
tegral has no speciality at all as compared to a general 
incommensurate filling. The Fermi surface has a almost 
perfect rounded shape. Thus, it is hard to guess what 
would happen in the strong correlation regime from weak 
coupling analysis. To the contrary, a quarter-filled sys- 
tem on triangular lattice with a negative hopping inte- 
gral is special in the sense that the chemical potential 
rests just on the Van Hove singularity of the density of 
state (/i = —2\t\ at zero temperature). The Fermi surface 
is nested with nesting vectors Qi,Q2, and Q3. As the 
result of the nesting property of the Fermi surface, the 
system develops magnetic order at an arbitrarily small 
U. In Ref. 12, the collinear state shown in Fig. la 
and the non-coplanar state shown in Fig.l.b are found to 
be the most favorable choice for the magnetic structure. 
Among the two, the non-coplanar state is found to be 
slightly more favorable. 

To have an idea on what would happen in the strong 
correlation regime for the quarter-filled Hubbard system 
with a positive hopping integral, we have conducted a 
unrestricted mean field search at zero temperature on a 
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12 X 12 lattice. The search is done in the mean field 
space with the local density rii and the three compo- 
nents of local spin density Si as variational parameters. 
We have used both the conjugate gradient method and 
the simulated annealing method to perform the search 
and used more than 100 random initial configurations 
for the local density and local spin density. We find the 
search always converge to the two configurations shown 
in Fig. 1 (up to global spin rotations and point group op- 
erations) for large enough value of y. Thus, although the 
quarter-filled system with positive hopping integral has 
no nesting property, it has the same magnetic ordering 
pattern in the strong correlation regime as the negative 
hopping system. 

The mean field equations for the collinear state is given 

by 



given by 



2U ^ 
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(Ck - E-),f{E-) 



{UmY 
, , (2) 

in which the sum is limited in the magnetic Brillouin zone 
and — — 2t(coska; + coskj^ + cos(ka. + kj,)) — /i is the 
dispersion on the triangular lattice. Here wc have used 
the convention for the momentum that k = k^jbi -|-k^b2. 
In this convention, the first Brillouin zone of the trian- 
gular lattice is given by (k^^'ky) G [— 7r,7r] ® [— 7r,7r] and 
the magnetic Brillouin zone of the collinear state is given 
by [0,7r] ® [-7r,7r]. E^ = g-+^-+Q3^Vfe-t.^+Q3)-+(^"-)- 
is quasiparticle energy in the magnetic ordered state. 

In the non-coplanar insulating state, the magnetic 

Brillouin zone reduces further into the region given by 
[0, tt] ® [0,7r]. The mean field equation for the non- 
coplanar state is given by 



(3) 



in which = ^Ea,^ «*(k)A^a,,3.^^(k) , a,^,n = 
1 • • • 8 and the sum is limited in the magnetic Brillouin 
zone. </'a(k) is the a — th component of the eigenvector 
for the mean field Hamiltonian with eigenvalue E]^. The 
mean field Hamiltonian matrix is of the following form 



H 



MF 



/ ai Ml 
Ml a+QiI 
MS 



M* 
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M2 

MS 
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M^ Ml Ck+Qal ) 



(4) 



here I is a 2 x 2 identity matrix, Mi 
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The matrix M. is 



M = - 
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f Ml M* M2 ^ 

Ml M2 M^ 

M^ M2 Ml 

\ M2 M^ Ml y 



(5) 



FIG. 2: The order parameters of magnetic ordered states at 
zero temperature as functions of the coupling strength U/\t\ 
for both signs of hopping integral. {a)t > 0, {h)t < 0. 

Fig. 2 shows the order parameter m as a function of 
U/\t\ at zero temperature for both signs of the hopping 
integral. For negative t , the magnetic order exist for 
all non-zero f as a result of the nesting property of the 
system. While for positive t, the magnetic order exists 
only for U/\t\ > Uj\t\ = 7.2. For both sign of t, the non- 
coplanar state is slightly more stable than the collinear 
, state at zero temperature. 

Fig. 3 shows the mean field phase diagram at finite tem- 
perature for both signs of the hopping integral. At fi- 
nite temperature, the gaplcss Fcrmion excitation in the 
collinear state will contribute more entropy than the 
gapped Fermion excitation in the non-coplanar insulating 
state. We find for sufficiently large U/\t\, the collinear 
state will become more stable than the non-coplanar 
state above a critical temperature at which a first order 
transition between the two occurs. In the phase diagram, 
the first order transition line between the collinear state 
and the non-coplanar state end at a critical end point on 
the magnetic-paramagnetic transition line. For positive 
t, the critical end point locates at U/\t\ = 7.9, while for 
negative t, it locates at a smaller value of U/\t\ = 5.5. 
Thus for both sign of the hopping integral, the collinear 
state is stable only in the strong correlation regime. 

It is interesting to note that the transition between the 
collinear state and non-coplanar state, though is of first 
order in nature, has almost no observable effect on the 
peak position and intensity of the magnetic Bragg peak. 
As we have explained above, the non-coplanar state has 
the same sets of magnetic Bragg peak as the collinear 
state, if the three domains of the latter as mixed with 
equal probability. In Fig. 4, we plot the change of the 
magnitude of the order parameter at the collinear - non- 
coplanar phase transition as a function of f7/|f|. We find 
the change to be always positive but its absolute value 
is always less than 1.5% of the ordered moment for all 
value of U/\t\ and for both signs of t. Thus we conclude 
that the neutron scattering experiments are not sensitive 
to the collinear-non-coplanar transition. 

However, as the system opens a full gap abruptly on 
the quasiparticle spectrum at the coUincar-non-coplanar 
transition, the longitudinal resistivity is expected to in- 
crease dramatically. Below the transition point, the sys- 
tem will enter the true insulating state. Thus in this sys- 
tem, the metal-insulator transition occurs in a two-step 
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manner. The quarter-filled system first breaks the lattice 
symmetry by a two-fold enlargement of the unit cell and 
enters a half-filled quasi-one dimensional metallic state at 
the paramagnetic-magnetic transition and then reduces 
its symmetry further by again a two-fold enlargement of 
the unit cell at the collinear-non-coplanar transition and 
enters a integer-filled insulating state. It is interesting 
to note that such a multi-step transition behavior is not 
limited to the quarter-filled system. At a more general 
filling fraction the system can first enter a ^(^)-filled 
intermediate metallic phase by a q{p)-tinies enlargement 
of unit cell and then arrive at the integer-filled true in- 
sulating state by a second transition with a p(g)-times 
enlargement of the unit cell. 

FIG. 3: The finite temperature phase diagram of the model 
for both signs of the hopping integral. (a)t > 0, (b)t < 0. 



FIG. 4: The jump of the order parameter at the collinear- 
non-coplanar transition point as a function of U/\t\ for both 
signs of the hopping integral. 

In the intermediate state between 53K and 88K, the 
collinear spin ordering induces dimensional reduction in 
the charge transport. In fact, all metal- insulator tran- 
sition can be viewed as a dimensional reduction process 
in the charge sector: an insulator can be viewed as a 
system in which the charge transport is quenched in all 
dimension. Such a dimensional reduction process in the 
charge sector can occur either directly, or like our exam- 
ple, in a multi-step manner in spatial dimension larger 
than 2. In the intermediate state , the charge trans- 
port is quenched in some dimension but remain active 
in other dimension. Such dimensional reduction in the 
bulk of the system can occur in condensed matter sys- 
tem with multiple low energy degree of freedoms. One 
famous example is the orbital ordering induced dimen- 
sional reduction in ruthenium oxides in which the elec- 
tron hopping is blocked except in one direction as a result 
of the symmetry property of the electron orbital wave 
function. Another example is the stripe-like structure in 
some cuprates superconductors in which the charge mo- 
tion becomes effectively one dimensional as a result of 
microscopic phase separation 14 1 . 

Finally, we discuss the relevance of our results to the 
physics of the NaQ,zGo02 system. The quasi-one di- 
mensional metallic state with collinear magnetic order 
has the proper characteristics to account for the exper- 



imental observations between 53K and 88K. Thus it is 
not unreasonable to assign this state to the intermedi- 
ate temperature region. However, as our model is purely 
two-dimensional, additional three dimensional coupling 
should be invoked to explain the finite temperature ex- 
istence of the magnetic order in the experiment. Here 
we assume that such three dimensional coupling will not 
alter the essential physics discussed in this paper. Fur- 
thermore, as our treatment of the model is limited to the 
mean field level, it is important to know how fluctuation 
effect will change the results. We leave this important is- 
sue to future study. As to the insulating state below 53K, 
we note that in real system, the perfect hexagonal sym- 
metry maybe broken by perturbation such as the sodium 
order potential, or spontaneously through the Peries in- 
stability process of the quasi-one dimensional intermedi- 
ate state, the real system have many different ways to 
become a true insulator. It is important to make sure 
that the non-coplanar state studied in this paper is sta- 
ble with respect to these perturbations before assigns it 
as the true ground state of the Nao,^Co02 system. 
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